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Abstract 

General aspects of the boundary value problem for the constraint equa- 
tions and their application to black holes are discussed. 



o 

1 Introduction 

o 

^J*- There are many kinds of boundaries in physics. For example, in electrody- 

namics, the boundary given by the interface between a charge distribution 
j^JT), and vacuum. However, this is not a fundamental boundary of Maxwell 

equations, in the following sense. This boundary is introduced in the 
sources by choosing a charge density with compact support. The sources 
, satisfy extra equations. Maxwell equations are expected to be fundamen- 

j.^ ■ tal equations; matter sources equations are phenomenological approxima- 

tion suitable to describe some specific matter models. The same apply to 
matter sources in Einstein equations. In this case the space time boundary 
is introduced in the sources of Einstein equations by choosing an energy 
density with compact support. 

In the case of vacuum Maxwell equations, the field does not interact 
with it self, hence it can not produce "its own boundary" . In the case of 
Einstein vacuum equation there exists such fundamental kind of boundary 
produced only by the self interaction of the vacuum field: black holes. 

In the context of an initial value formulation, the first step in order to 
understand the space time boundary produced by a black hole is the study 
of the intersection of this boundary with a spacelike three-dimensional 
Cauchy hypersurface. That is, to study the black hole boundary value 
problem for the constraint equations. This problem has been recently 
studied in [TTJ 

Since only fundamental properties of gravity are involved in a vacuum 
black hole, it can be expected that black hole boundary conditions can be 



1 



S. Dain 



2 



written in a geometrical form. It turns out that this is true. Moreover, 
black holes boundaries for the constraint equations suggest a deep inter- 
play between Riemannian geometry, elliptic equations and physics. In the 
present article we discuss some general aspects about this interplay. 

In section |3 we present the constraint equations and the correspond- 
ing boundary conditions. In section|2]we discuss elliptic reductions to this 
equation. General kind of boundary conditions are discussed in section [I] 
In section |5] we discuss boundary conditions that are physically meaning- 
ful. Finally, in section [|3 we discuss black hole boundary conditions. 

2 The constraint equations 

Let Ck be a finite collection of compact sets in R 3 . We define the exte- 
rior region tl = R 3 \ U^Cfe. For simplicity, we will mainly consider the 
constraint equations in the time-symmetric case, that is, when extrinsic 
curvature vanishes. However, all the following consideration apply also 
to the general case (see Let h a b be a Riemannian metric on Cl and 
let R be the corresponding Ricci scalar. The time-symmetric, vacuum, 
constraint equation is given by 

R = 0, (1) 

on fi. 

There exists two kinds of boundary conditions for equation Q in fi: 
outer boundary conditions and inner boundary conditions. The outer 
boundary condition is asymptotic flatness, and it is essentially a fall of 
condition on h a t- Physically, it means that we have an isolated system. 
The data will be called asymptotically flat if there exists some compact 
set C, with UfcCfe C C, such that Q\C can be mapped by a coordinate 
system x? diffeomorphically onto the complement of a closed ball in R 3 
and we have in these coordinates 

h ii = {l+-^)S iJ +0(f-% (2) 

as f = (X^=i (x-') 2 ) 1 ^ 2 — > oo, where the constant m is the total mass of 
the initial data. 

The inner boundary condition will be the black hole boundary con- 
dition. The boundaries dCk are assumed to be smooth, two dimensional 
surfaces in (fi, h). Let u a be the unit normal of dCk, with respect to h a b, 
pointing in the outward direction of Cl. Let t a be the unit timelike vec- 
tor field orthogonal to the hypersurface £1 with respect to the spacetime 
metric g a b (t a t b g a b = — 1 with our signature convention) The outgoing 
and ingoing null geodesies orthogonal to dCk are given by l a = t a — v a 
and k a — t a + v a respectively, the corresponding expansions are given by 
Q+ = V a l a and 0_ = X7 a k a , where V a is the connexion with respect to 
g a b- We can calculate these expansions in terms of quantities intrinsic to 
the initial data. In the particular case of time-symmetric data we have 



e_ = -e+ = ii, 



(3) 
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where H = D a i> a , and D a is the covariant derivative with respect h a b- 
That is, H is the mean curvature of the two-dimensional surface <9f2 with 
respect to the metric h a b and the normal v a . The boundary will be called 
marginally trapped if 

H = 0. (4) 

A marginally trapped surface indicates the presence of a black hole (see 
|14) and also the discussion in 9 ) . Equation JJJ geometrically means that 
the two dimensional boundary is an extremal surface with respect to the 
metric h a b- Equation gj will be our inner boundary condition. 

3 The constraint equations as an elliptic 
system 

We want to find solutions h a b to equation (QQl which satisfy boundary 
conditions J3J and Q. If we write equation Q in terms of the metric 
components of h a b we get a complicated non linear equation. We have six 
unknown functions in the metric h a b and only one equation. That is, we 
have an underdetermined system. The strategy to solve this equation is to 
split the six unknowns into two sets. One set will be called the free data, 
we want to prescribe them freely or at least with some restrictions easy to 
achieve, which, in particular, do not involve solving differential equations. 
Roughly speaking, the free data set should contain five free functions. The 
other set will contain only one function. For a given choice of free data 
we have to solve equation (0 to calculate this function. 

There exists a priori many ways of doing this splitting. For example 
we can chose some coordinate system and chose the free data set to be 
five components of the metric in these coordinates and try to solve the 
equation for the remainder. But the resulting equation will be in general 
a complicated non linear equation which is of no known type and hence 
will be difficult, if not impossible, to prove that in fact for every choice of 
free data we do get a solution for the remainder. 

In order to control the behavior of the solution h a b at the boundary we 
need also to prescribe boundary conditions for the unknown function. Our 
splitting and boundary conditions will be successful if for arbitrary free 
data and boundary conditions we always get a solution for the reminder. 
This suggests that the problem has an elliptic nature. To some extend 
this is true, the constraint equation can be reduced to an elliptic system. 
However, this is not the only way of solving them, for example they can 
be reduced to a parabolic system [3] and this lead to the discovery of 
new kinds of solutions. Nevertheless, so far only the elliptic approach has 
been successful in the general case (the parabolic system has been mainly 
studied in the time symmetric case). 

There is not a unique way of getting an elliptic system out of equation 
Different elliptic systems will lead to different choices of free data 
and boundary conditions, and hence they can be more appropriate to 
describe different kinds of physical situations. Here we will use the so 
called conformal method, which is probably the simplest one, since it lead 
to semilinear equations, which reduce to a linear equation in the time- 
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symmetric case. Other elliptic reductions give quasi linear equations [7| 
OH El 011]' So far, black hole boundary conditions have been only studied 
using the conformal method (11) 0. It will very be interesting to study 
them with other elliptic reductions. 

The conformal method is as follows. Let h a b be a Riemannian metric, 
let t/)bea positive solution of the following equation 



where Lh = D a D a — R/8, D a is the covariant derivative with respect to 
hab and R is the Ricci scalar of the metric h a b- Then, the rescaled metric 
h a b = ip^hab will satisfy equation Q. Note that the differential operator 
Lh is elliptic, hence the linear equation Q is an elliptic reduction to the 
time symmetric constraint equation Q. 

Two metrics h a b and h a b belong to the same conformal class if there 
exists a positive conformal factor ip such that h a b = tjj 4 h a b- For any 
metric in the same conformal class we get the same solution h a b- That is, 
the free data set is given by the conformal class of metrics, which can be 
represented by five free functions. 



What kind of boundary conditions are compatible with the constraint 
equations? If we have reduced them to an elliptic system, this question 
can be answer in full generality. Given an elliptic system, the boundary 
conditions will lead to a well posed problem if and only if they satisfy 
the so called complementing condition or Lopatinski-Schapiro conditions 
(see pQ |5] for a precise statement and also the introductory book |13p. 
These are conditions at the linear level. For non linear system, these 
conditions are imposed to the associated linearized problem. Of course 
in the non linear case, these conditions are in general only necessary but 
not sufficient to prove the existence of solution. For non linear systems we 
have to study each particular case to decide whether there exists solutions 
or not. 

Simple examples of boundary conditions that satisfy the Lopatinski- 
Schapiro requirements for equation Q are Dirichlet and Neumann bound- 
ary conditions for the conformal factor ij). More general boundary condi- 
tions are possible, it is even possible that the order of the derivatives in 
the boundary operator is higher than the order of the derivatives in the 
differential operator. 

However, here we are interested in positive solutions. This is an extra 
requirement of our particular elliptic reduction. If ip is zero at some point 
then hab = ^h a b will not be a Riemannian metric at this point. The 
positivity of the solution can be proved by a particular feature of second 
order elliptic equations: the maximum principle (see, for example, |1U|). 

In order to use the maximum principle we need an extra requirement 
on the lower order coefficients of the operator Lh 



L h ^ = 0, 



(5) 



4 Elliptic boundary conditions 



R > 0. 



(6) 
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The most general kind of boundary condition that satisfy the Lopatinski- 
Shapiro conditions and also allow us to use the maximum principle to 
prove positivity is given by the Dirichlet boundary condition 

ip = p> on <9fi (7) 

or the oblique derivative boundary condition 

p a D a il> + cap = tp on 9fi, (8) 

where p, a and /3 a are arbitrary functions which satisfy p > 0, a > and 
P a u a > on dQ., where u a is the outward unit normal with respect to the 
conformal metric h a t- We note that p a v a > guarantee that JSJ satisfies 
the Lopatinski-Schapiro conditions. 

In the exterior region Q we need in addition the asymptotic flatness 
condition which in this case is given by 

lim if) = 1. (9) 

r — >oo 

Conditions JSJ, Q or JBJl, @ will guarantee the existence of a unique 
positive solution ip of equation and hence an asymptotically flat metric 
h a b which satisfies the constraint equation Q with or © as inner 
boundary condition. Of course, for arbitrary p, a and f5 a , the boundary 
condition (JHJ will not have any interesting geometrical meaning. 



5 On physical boundary conditions 

Probably any smooth solution of the constraint equation in some bounded 
region can have some physical interpretation in the sense that this bounded 
region can be a piece of a space time that can describe some physical phe- 
nomena. The situation is different in the case of an exterior region f2: in 
order to be physically meaningful as a description of an isolated system 
the solution must have positive mass. Since we artificially cut out a region 
in R 3 , the positivity mass theorem does not automatically apply. Many of 
the solutions found in the previous section will have negative total mass. 
This can be explicitly seen in the following example. 

Let us consider the Schwarzschild, time symmetric, initial data. In 
this case h a ^ = S a ^, where 8 a ^ is the flat metric, — A, where A is the 
flat Laplacian and the conformal factor is given by 

tf = l + £. (10) 

We chose the exterior region Cl to be the exterior of a ball of radius r = a. 
Note that with our conventions v a = ~(d/dr) a . The Dirichlet boundary 
condition is given by 

ij) = ipo on dCl, (11) 

where po is a positive constant. Using (I10P one easily check that po < 1 
implies m < and po > 1 implies m > 0. Take a > ra/2, then we will 
have ip > in the exterior region and m < if we chose ipo < 1. This 
means that p>o < 1 is not a physical boundary condition although it math- 
ematically consistent in the sense that it gives us existence and uniqueness 
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of a positive solution. In the general case, if we impose Dirichlet boundary 
conditions to the conformal factor appears to be difficult to recognize for 
which value of the boundary function we will get data with positive mass. 
We conclude that Dirichlet boundary conditions for ip are not physically 
meaningful in general. 

In this example, the Neumann condition 



a n , dip 
or 



= <fio (12) 



with tpo > always produce data with positive mass, since m = 2(poa 2 . 
This can be generalized. Consider the case of non trivial extrinsic curva- 
ture, and assume that the trace of it is zero. Assume that the conformal 
metric h a b satisfies R = and the following fall off hij = 5ij + 0(r~ 2 ). 
The equation for the conformal factor is given by 

L h ^ = D a D a i = -^^, (13) 

where K a b is the rescaled extrinsic curvature. Let us impose Neumann 
boundary condition to this equation 

u a D a i = tp on an, (14) 

with ip > 0. Integrating equation 11311 on Cl and using the Gauss theorem 
we get 

m> I u a D a i>dS > 0, (15) 
where we have used 1141 and the following expression for the mass 

m=--L lim / ^-dS r . (16) 

4-7T i — >oo J s or 

We conclude that every data which satisfy the boundary condition 114H 
in the appropriate conformal class will have positive mass. Note that the 
boundary condition 1141 is linear even in the general, non time-symmetric, 
case. It is not clear to me the meaning of this condition, and I am not 
aware of any application of it. Although in |11| only black hole boundary 
conditions has been studied, remarkably, all the solutions founded there 
satisfy in addition condition ffll4B . 



6 Black hole boundary conditions 

The most important inner boundary condition for the vacuum constraint 
equation is the black hole condition This boundary condition has 

both a geometric and physical meaning. There exist versions of the pos- 
itivity mass theorem which include black holes as inner boundaries (see 
for example |12p. 

Condition Q can be written in terms of the conformal quantities as 



4v a D a i/j + H^ = 0, 



(17) 
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where H = D a u a . 

If H > on <9f2, then condition l|17|l has the form (Hj. If, in addition, 
the conformal metric h a b satisfies R > on f2 we can apply the standard 
linear elliptic theory and the maximum principle to prove that there will 
exist a unique positive solution (j> of our problem. 

The conformal metric is not really a free data, it should satisfy the 
condition R > 0, and also the boundary dQ should satisfy H > 0. Note 
that if we chose the conformal metric to be the flat metric and dfl any 
sphere, then this boundary will have H < 0, and hence it does not sat- 
isfy our hypothesis. However, it is simple to construct families of metric 
which satisfy R > and H > on dCl. Consider, for example, the time- 
symmetric initial data for Reissner-Nordstrom. The metric is given by 

hab = ffiSab, (18) 

where 

i> = — \/{q + 2r + m)(-q + 2r + m). (19) 
Ir 

The constants q and m are the charge and the mass of the data respec- 
tively. We assume q 2 < m 2 . The Ricci scalar is given by 

(20) 



which is positive. For r < ro — {\/m 2 — Q 2 )/2 the mean curvature of the 
two surfaces of constant radius is positive. Then this metric satisfies our 
hypothesis. Moreover, take q > (this implies R > 0), let dfl such that 
H > 0, let e small enough and X a b and arbitrary tensor field; then the 
metric h ao + e\ a t satisfies also R > and 9f2 will satisfy H > with 
respect to this metric. 

In the general case, black hole boundary conditions are non linear. 
This introduces extra difficulties in both the existence proof and the elec- 
tion of free data. This problem has been recently studied in II 1 1 In 
those references, large classes of black hole exterior regions have been 
constructed. However, it is still an open problem how to construct and 
characterize all possible initial data for black holes exterior regions. 
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